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Presents a mathematical jet engine where 
the measured transition functions of the engine thrust can 
be simulated in a quasi-stationary operation as well as in 
acceleration or deceleration schedules. Because of its sim- 
plicity and high fidelity, it is especially suited for rep- 
resenting jet engines in digital simulation programs. 
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Short Communications 

A SIMPLE PROPULSION SYSTEM MODEL FOR THE SIMULATION OF 
NONLINEAR DYNAMIC THRUST RESPONSE. 

■ 
Gunther Schanzer and Peter Krauspe, TU Braunschweig 
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Summary: A simple mathematical jet engine description is presented 
with which the measured transition functions of the engine thrust 
can be simulated in quasi-stationary operation as well as in 
acceleration or deceleration schedules. Due to its simplicity 
along with high fidelity, it is especially suited for. the 
representation of jet engines in digital simulation programs. 



1 . Introduction / 

The representation of measured thrust-time sequences of jet 
engines in digital simulation programs, of whose accuracy particularly 
high demands are made whenever the jet engine thrust is the primary 
control quantity of a regulation process, often presents problems. 
For one, the transition functions (figure 1) can be represented in 
approximation only by linear control paths of very high order 
( 2C and more) on account of the large ratios of delay time and 
recovery time. For another, the transition behavior of the thrust 
changes non-linearly as a function of the engine rpm (figure 2). 
Finally the time behavior, because of internal energy limitations 
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of the allowable air mass flux, is dependent on the magnitude of 
the control quantity jump; here one must differentiate betv/een 
operation on the so-called steady-state operating line (quasi- 
stationary operation) and the so-called acceleration/deceleration 
schedule. The adaptation of a linear control path of high order 
to the named condition is very costly Cl]« The description of the 
engine dynamics in the phase plane £2] here offers the possibility 
for a simple and clear representation of the linear and non-linear 
portion. The model parameters used depend on engine- and flight 
conditions, such as pressure ratios, air mass throughput, inlet 
temperature, Mach number, etc. For reasons of simplicity, 
one must knowingly forego to include the physical relations between 
the engine parameters into the model, as is generally customary 
in the simulation of the entire engine C3l- The model parameters 
of the presented model are gained basically from thrust measurements 
(ground station test or simulation), which apply for the particular 
flight condition. 

2. Symbols 

F thrust 

p thrust change with time (= dF/dt) 

F-r idle thrust 

F „ maximum thrust 

t time 

T. time constant 

Ng rpm of the high pressure compressor 

ftp rpm changes with time (= dNp/dt) 
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Figure 1: Transition behavior of 
different jet engines. 
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Figure 2: Jump responses of a 
jet engine (measurements) for 
a) large and b) small thrust 
changes. 
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Figure 3 s Phase diagram of the 
thrust sequence of a jet engine 
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Figure 4-: Comparison of measured 
and simulated thrust curves . 



3. Bases and physical relations 

The dynamic non-linear engine "behavior, such- as is presented, 
e.g., in figure 2, can be approximated (figure 3) with good accuracy 
by straight lines at least for the multi-shaft engines investigated 
JT3D, JT9D, CF6-50C in the phase plane P, F. The desired function 
F(t) can, since an F value is assigned to every value of F, be gained 
from a simple integration of F: 

F{t) = F + JF(F)dt. (1) 





In figure k a plot of the thrust-time curves calculated in this / 
way is plotted along with the measurements on which they are based. 
As the result of the additional statement of the rpm N~ of the high 
pressure compressor, the typical dynamic properties of a jet engine 
power plant become more pronounced (figure 5) • If one subdivides 
the phase curves in figure 3 ssid the assigned jump response to a 
change in gas throttling in figure 5 into subregions, one can find 
a clear physical interpretation for every individual region: 

(I) Starting region 

Increase of the engine speed N ? from idle with nearly constant, 

relatively small changes in the speed #2.,. • One must remember 

that for turbine- jet—stream engines there is a strong non-linear 

relation between the thrust and the speed. Thus, one must assign, 

'e.g., to the idle thrust F T an idle speed of nearly 60% of N, : 

jj c ,max 

thus, the entire thrust spectrum from F T to F lies within the 
. ii max 

remaining ^-O^.of the allowable speed range. In region I one obtains, 
despite the non-linear behavior of the speed, a nearly constant 
value of F. • 



K 



(II) Transition region to maximum acceleration N Q mnv 

c. i max 

In the phase plane the slope F/F of the phase curve is constant. 

(III) Region of constant thrust change F with almost linear 
dependency on thrust and speed 

The acceleration of a multi-shaft jet engine is limited in 
this range by the pump limitation of the high pressure compressor. 
In multi-shaft engines the safe distance from the pump limit is 
guaranteed by engine controllers. Thus, the dynamic behavior and 
particularly the maximum acceleration capability is determined in 
region III by the combined action of engine and controller. 

(IV) Transition region to the stationary final value 

The time behavior corresponds to the behavior of a linear 
system of first order with the characteristic time constant T 
(figure 6) : 

'f(0 = ^Nom.W-F<0]=^r-AF(0. (2) 

The slope of the phase curve in the phase plane thus is given by 

• (^),v=T =C ° nSt - (3) 

Let the characteristic time constant of region IV be T ? . If one 
compares the transition regions II and IV one can find a formal 
agreement: both sections were given in the phase plane by a straight 
line slope F/F. Let the characteristic time constant T-, be defined 
in region II analogous to equation {"})'• 



(-) = -- 

\AF/n T, 



const. 



(1) 
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Figure 5* Relations of thrust and speed of a multi-shaft 
jet engine . 



A comparison with figure 5 shows that for the jet engine investigated 
here one can use the approximation 

T T T . (5) 

These time constants characterize time behavior of the jet engine 
for small characteristic value jumps ("small time behavior") and 
lie in the range from 0.25 to 1 second. However, one must note 
that the model parameters F f T,, , T n and T 9 apply each time only 

IT13.X V X *C 

for a given ambient- and flight condition. The model parameters 
presented here refer to ground stand tests under normal atmospheric 
conditions. As is known from flight measurements, larger values 
for T, and Tp are obtained at higher altitudes as the result of 
decreasing air density and the reductions in air mass flow resulting 
therefrom. 



With the aid of the small-time behavior one can also describe 
i 
those steady state changes which are characterized by relatively 

small control quantity jumps, which the system can follow without 

large delays. The characteristic time constant T of the linear 

system can be obtained by tangent formation directly from the 




Figure 6: Jump response and phase diagram of a linear system 
of first order. 



measured curve of the transition function of the path, and namely, 
as shown in figure 2, from the jump response to a) large as well 
as b) small jumps: 



T a ■ T b = T v 



06) 



It is essential that the path can take on practically all possible 
intermediate values of F v/hich lie below the maximum value F 
of the acceleration line. Thus, the steady-state behavior of the 
engine is obtained in the phase plane by the number of all points 
(F, F) within the phase curve. 
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Figure "?: Block switch diagram of the engine model. 
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Figure 8-: Examples of thrust commands in the phase' plane. 

a) 75 Wo (steady-state) ' 

b) 40 50$ (limited acceleration) 



3. Model description 

The described relations are shown in the block switch diagram 
in figure 7* According to equation (2). one can determine as a 
function of the difference between desired- and actual value of 
ths thrust F the initial value F of a linear system with the time 
constant T . If the thus determined point (F, F) lies in the phase 
plane within the surface bounded by the phase curve (figure 8a), 

• 

one is dealing with a change of state, and the desired value F 
is equal to the result from equation (2); if the point lies outside 
of the phase boundary (figure 8b), one must put in for F the 
maximum possible thrust change corresponding to the associated 
value of the ordinate of the phase curve. In figure 7 this action 
is achieved by the function of the limiting device B whose charac- 
teristic line is the phase curve of the jet engine. A subsequent 
integration of F in accordance with equation (1) provides the 
desired value,F(t). The same is true for the operation of the 
deceleration line or in the region of the steady-state deceleration 
of the jet engine. 
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The good agreement of the measured and simulated thrust 
curves, together with the low programming expenditure, seemed 
to make, the presented jet engine model to he very useful. 
Applications with detailed numerical calculations are contained 
in M . 
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